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The dipole P{F) of systems with periodic boundary conditions (PBC) in a static electric field F 
is applied to one-dimensional Peierls-Hubbard models for organic charge-transfer (CT) salts. Exact 
results for P{F) are obtained for finite systems of A'^= 14 and 16 sites that are almost converged 
to infinite chains in deformable lattices subject to a Peierls transition. The electronic polarizability 
per site, a^i = {dP/dF)o, of rigid stacks with alternating transfer integrals t{l it S) diverges at the 
neutral- ionic transition for 5 = but remains finite for 5 > in dimerized chains. The Peierls or 
dimerization mode couples to charge fluctuations along the stack and results in large vibrational 
contributions, a^it, that are related to dP/dS and that peak sharply at the Peierls transition. 
The extension of P{F) to correlated electronic states yields the dielectric response k of models 
with neutral-ionic or Peierls transitions, where k peaks > 100 are found with parameters used 
previously for variable ionicity p and vibrational spectra of CT salts. The calculated k accounts for 
the dielectric response of CT salts based on substituted TTFs (tetrathiafulvalene) and substituted 
CAs (chloranil). The role of lattice stiffness appears clearly in models: soft systems have a Peierls 
instability at small p and continuous crossover to large p, while stiff stacks such as TTF-CA have a 
first-order transition with discontinuous p that is both a neutral-ionic and Peierls transition. The 
transitions are associated with tuning the electronic ground state of insulators via temperature or 
pressure in experiments, or via model parameters in calculations. 



I. INTRODUCTION 

Torrance et ali reported the first neutral-ionic transi- 
tion (NIT) in the organic charge-transfer (CT) salt TTF- 
CA. Tetrathiafulvalene (TTF) is a potent 7r-electron 
donor (D), while chloranil (CA) is a strong 7r-acceptor 
(A). As confirmed subsequently by detailed vibrational 
and structural studiesf?*^ the 300 K structure contains 
mixed regular stacks ..D''+A''^D''+A''~... with face-to- 
face overlap and ionicity p ~ 0.3. Equal intermolecu- 
lar interaction with each neighbor follows from inversion 
symmetry at the centers of TTF and CA. Below the first- 
order transition at T^ = 81 K, the dimerized ionic stacks 
have /5 ~ 0.7 and alternately short and long intermolec- 
ular spacing. The ionic lattice is unstable to a transition 
that becomes a spin-Peierls transition^ in the limit p —> 
1. In TTF-CA and other CT salts with discontinuous p 
at Tc, the neutral-ionic and Peierls transitions coincide. 

Systems with neutral-ionic or valence transitions are 
necessarily rare: two ground states (GS) with differ- 
ent charge densities must be almost degenerate in or- 
der to switch between them by changing temperature or 
pressure. Horiuchi et aliSi^iSiiSiii have recently stud- 
ied NITs in a series of substituted TTFs and CAs whose 
static dielectric responses have large peaks (k > 100). 
Their interpretation is a quantum phase transition be- 
tween a neutral-regular and an ionic-dimerized GS con- 
trolled by an external parameter such as temperature or 
pressure,.-^ ^ No modeling of k is proposed, however, and 
no modeling is possible without the static polarizability 
of extended systems. In this paper we compute the static 



polarizability of Peierls-Hubbard models for organic CT 
salts and vary GS parameters to induce transitions. We 
obtain explicitly electronic and vibrational contributions 
to K for correlated states. Systems with continuous p 
at the K peak and p < 1/2 throughout are of particu- 
lar interest, since they do not fit the simple paradigm of 
neutral-regular or ionic-dimerized. 

The dielectric constant, k, is the GS polarizability a, or 
induced dipole moment, per unit cell in the crystal. The 
modern theory of polarization in insulatorsi^ has resolved 
the problem of a position or dipole operator in systems 
with periodic boundary conditions (PBC). Resta^'^ intro- 
duced the following expression for the expectation value 
of the dipole moment, P: 

P= — Jm ln(^|exp(i )\tjj) = — Im\nZ (1) 

27T N 27r 

Here is the exact electronic GS, M is the conven- 
tional dipole operator and PBC are ensured via super- 
cells. Polarization in insulators is related to the phase 
rather than the amplitude of and can be formulated as 
a Berry phase. P applies to correlated as well as non- 
interacting systems and provides a unified approach to 
previous treatments of polarization and GS localization 
in insulators.^'* The magnitude of Z is finite in insulators 
and zero in conductors. The phase or twist operator in 
Eq. (1) also appears in other contextSfi^ We will work 
with P for the exact GS of modified Hubbard models 
with N sites. 

Although a static electric field F is not compatible 
with PBC, Nunes and Gonze^^ showed recently how to 
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combine F and P. They propose to define P{F) with 
')Jj{F) in Eq. and to minimize a functional that, in a 
one-dimensional system, is 

E{i;F, F) = EO^f) - NFP{^f) (2) 

Several groupsiiiSiiSiSS have applied Eq. Q to the elec- 
tronic dielectric constants of metal oxides and silicates, 
using Wannier functions for 'il'{F) and perturbation the- 
ory to find corrections to "0(0) • As shown below, cor- 
related models with finite N are simpler because the 
basis is finite, which allows exact solution for V'(-f') iii 
Eq. The procedure in Section 2 differs from the 

model-exact treatment of static or dynamic nonlinear op- 
tical responseSfSi however, since the dipole operator M 
is now in the exponent. 

The optical phonon of the mixed regular stacks 
..'DP'^ AP^DP'^ ... describes a rigid displacement of ei- 
ther sublatticc, and in the hypothesis of linear elcctron- 
phonon (e-ph) coupling leads to alternating CT integrals: 
t(l+(5), t(l— (5). We assign the optical phonon a harmonic 
potential with frequency at p = 0. With increasing 
p, e-ph coupling generates large anharmonicitiesS^ and 
an overall softening of the vibrational frequency, which 
vanishes at the Peierls transition. A static electric field, 
F, strongly affects the lattice degrees of freedom at in- 
termediate p, especially near the Peierls transition. The 
polarizability per site has then substantial electronic and 
vibrational contributions, 

dP dP dP 85 

The P derivatives above illustrate that physical processes 
are associated with changesi^ of P and underscore the 
central role of P for the dielectric response. In the vibra- 
tional part, ^ is related to the IR intensity of the Peierls 
moda^ and ^ to the field-induced change in dimeriza- 
tion. As anticipated from non-interacting systemsf22, ami, 
has a peak at the Peierls transition that is further en- 
hanced by the frequency softening. The electronic part, 
by contrast, diverges at the NIT of the rigid regular lat- 
tice but varies slowly near the Peierls instability of de- 
formable lattices. 

The paper is organized as follows. Section 2 presents 
the exact solution of P{F) using Eq. 10) for quantum 
cell models with a large but finite basis. The electronic 
polarizability of modified Hubbard models with fixed al- 
ternation is obtained in Section 3. The roles of lattice 
vibrations are introduced in Section 4 and developed ex- 
plicitly for the dimerization mode and its softening at 
the Peierls transition. Section 5 contains polarizability, 
ionicity and related results for Peierls-Hubbard models 
with equilibrium dimerization. The k peak at the Peierls 
transition is due to amb and increases with the lattice 
stiffness. We discuss in Section 6 the applicability of 
model results to dielectric peaks in organic CT salts. 



II. STATIC ELECTRONIC POLARIZABILITIES 

Quantum cell models typically have a single Wannier 
orbital (pp per site that, in principle, can incorporate 
intra-site correlations.^^ On-site repulsion [/ > is kept 
in Hubbard models and long-range Coulomb interactions 
in the Pariser-Parr-Pople model. Total spin S is usually 
conserved, as in models for CT salts, but this condition 
can be relaxed. The defining feature of cell models is a 
large but finite many-electron basis in real space that is 
complete for a finite number of orbitals 4>p. Slater deter- 
minants of 0p's can be used in general for eigenstates of 
5"^, but linear combinations of determinants are required 
for the smaller basis when S is conserved. Easily visual- 
ized valence bond (VB)2^ diagrams are then convenient 
and used below. 

The normalized singlet [S — 0) diagram \k^q) speci- 
fies the occupation number = 0,1 or 2 at each site 
p = 1,2, ...TV and the spin pairing of all singly-occupied 
site. The index q is introduced below for the eigenvalue 
Mk of the dipole operator. The VB basis is orthogonal 
with respect to charge, but there is overlap Suk' between 
diagrams with identical Upk and different spin pairing. 
The operator M for N sites with uniform spacing a = 1 
and e = 1 is 

N N 

M ^^PQp P{zp - Up) (4) 

p—i p—i 

Here qp = Zp — Up is the charge operator, with Zp = at A 
and Zp = 2 at D, and there are two electrons per DA unit 
cell. The neutral GS has regular spacing a along the stack 
while the dimerized ionic GS has alternating spacing 
a±u. For simplicity we use M in Eq. Q for both regular 
and dimerized GS, thereby neglecting the motion of ±ep 
charges; corrections for ±u can readily be includedj22, M 
is diagonal in the \k^q) basis. Its eigenvalues are the in- 
tegers, modulo N, Mk = 0, ±1, ±2, .., -iV/2. We take 
them as q = 2TTMk/N for — tt < g < tt in the first Bril- 
louin zone. The matrix element in Eq. is thenSS 

(V'lexpi^^^l^) = ^ W,expig= Zj. + Zj, (5) 

where Wq > is the weight of basis vectors \k,q) with 
q — 2TTMk/N . Diagrams with q ^ Q have greatest weight 
on the neutral side and lead to Z^. ~ 1, while diagrams 
with q ^ ±n and —1 dominate on the ionic side. 

The regular chain has real Z by symmetry and Z — at 
the NIT. 

Since the VB basis is complete, the GS can be ex- 
panded as 

^(i^)=^Cfe,(F)|fc,g) (6) 

k,q 

The linear coefficient Ckq are real and known at = 
by hypothesis. Substitution of Eq. (jSjl into Eq. Q and 
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varying the Ckq leads to a secular determinant. As dis- 
cussed elsewherCfS^ sparse-matrix methods that exploit 
completeness are much more convenient for large ( > 10^) 
basis sets. We consequently seek an operator whose ma- 
trix elements correspond to P{F) in Eq. The partial 
derivative of P{F) with respect to Ck'q is 



dP{F) 



'q 



k.q 



Zx (F) sin q — Zy (F) cos q 
2n[Z.AFr + Zy{Fr] 



(7) 

We define the quantity in {...} as /ifc^ and interpret it as 
the induced dipole of diagram \k,q). Since precisely the 
same result is obtained by introducing diagonal energies 
—NF/jLkq for each diagram, the desired induced-dipole 
operator is 



^ ' 1-K Z{F) 



(8) 



It follows immediately that \k,q) is an eigenfunction of 
AM with eigenvalue Nfikq, as given by Eq. ijJJ). The large 
response or polarizability for small Z{F) is also evident. 

The minimization of the energy functional in Eq. (j^J at 
finite F is equivalent to finding the GS of the Hamiltonian 



H{F) - FAM{F) 



(9) 



Given the action of H{F) on any diagram, we obtain 
ip{F) as usual^^ except for one complication. Since we 
know Z(0) from Eq. |(TJ, but need Z{F) in Eq. ©, it- 
erations are needed. We start with AM(0) to obtain 
Z^-^^F) and hence AM^^^F), which is then used in 
Eq. © to find Z^^\F) and AM(2)(F). We repeat un- 
til Z{F), iP{F) and P{F) have converged. The GS with 
AM(0) requires the same computational effort as V'(O) 
and suffices for the linear response, aei = -P'(O). More 
precisely, we use AM(0) in Eq. © to find ip{F) at small 
F, compute P{F) in Eq. Q and evaluate P'(0). We 
have {tplF)\AM{0)\ip{F)) = FP'{0) and, as seen directly 
from Eq. iQ, {2lj{F)\AM{F)\il;{F)) = The interpre- 
tation of Eq. © as the induced-dipole operator is now 
clear: AM(0) generates the linear correction to P{F). 

We note that the Berry-phase formulation of P{F) in 
Eq. ^ differs from standard expressions for finite sys- 
tems with open boundary conditions. In finite systems, 
the interaction with a field is H — FM instead of the 
PBC Hamiltonian in Eq. The dipole operator M is 
diagonal in the VB basis, with — J2pP1pk for dia- 
gram \k,q). Now ^{F) is the GS of iJ - FM and the 
electronic polarization is ($(F)|M|$(F)). In contrast 
to PBC, there is no iteration and Mk is not modulo N. 
The theoretical ideas and computational difhculties that 
motivated the PBC formulation of P in Eq. Q and its 
coupling to F in Eq. ^ do not arise in quantum cell 
models, whose limited basis is unphysical in this respect. 
But extrapolations of linear and nonlinear optical coef- 
ficients to infinite systems have been widely discussed~ 
for both noninteracting (Hiickel) and correlated (Hub- 
bard, Pariser-Parr-Pople) models, especially in connec- 
tion with frequency-dependent responses. The velocity 



operator version of M is often used to estimate transition 
dipole moments for extended (PBC) systems and is con- 
sistent with the real-space resultsiSl Susceptibilities ob- 
tained with the familiar sum-over-states (SOS) approach 
require all excited states, a difficult task for systems with 
large correlated basis, but can be related to GS responses 
with correction vectors. The calculation of P{F) and its 
derivatives only requires the GS. 



III. ELECTRONIC POLARIZABILITY OF 
MODIFIED HUBBARD MODELS 

To illustrate the PBC polarizability of rigid lattices, 
we consider Hubbard models that have been applied to 
the electronic structure of both organic CT salts^^'^^ and 
transition-metal oxidea^ii. Face-to-face 7r-overlap in CT 
salts leads to transfer integrals t between sites with en- 
ergy — A at D and A at A. The modified Hubbard model 



-28 



^Upijlp - 1) 



E 



A{~lYn.p + U- 



(10) 



We take i = 1 as the unit of energy, consider alternation 
< (S < 1 along the chain, and assume equal on-site 
U > Q for A and D sites. Coulomb or other interactions 
can be added to Hhu without increasing the correlated 
basis. Such models conserve S, have Cjv/2 symmetry and 
require comparable effort to obtain the exact GS. 

On physical grounds, high-energy and A^~ sites 
can be neglected to obtain a restricted basis that is use- 
fully smaller. Formally, we define F = A — C//2 and take 
the limit A, C/ ^ cx) in Eq. H10|) while keeping finite F. 
The electronic problem for rigid lattices simplifies to 

p.cr 

+ Y^T{-irnp (11) 

p 

The GS for F >> is the neutral (p 0) lattice with N 
electrons paired on N/2 donors, while the GS for F << 
is the ionic lattice (p ~ 1) with one electron per site and 
overall singlet pairing. The regular {S — 0) stack has an 
NIT as a function of A in Eq. ^ or F in Eq. ^ that 
has been studied by different methods, including exact 
GS energies and p up to = 16 for Eq. (|10|l and up to 
N ^22 for Eq. 

In the present work, we take ^'(0) in Eq. JHJl to be 
the exact GS of HoiS,T). P(0) in Eq. (jTJ is the relevant 
GS dipole per site in units of ea. Then, using AM(0) 
in Eq. JSJ, we obtain ^{F) exactly, compute Z{F) and 
P{F), and evaluate a^i = ^'(0) numerically. Figure^ 
and b show, respectively, aei and Z(F) for 6 = stacks of 
A^ = 14 (circles) and 16 (crosses). To understand these 
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FIG. 1: The logarithm of the hnear electronic susceptibility, 
Qfe! , and the modulus of Z vsF calculated for the Hamiltonian 
in Eq. Illlll with 5 = for stacks of 14 and 16 sites (circles 
and crosses, respectively). The arrows mark the NIT of the 
infinite strand. 



results, we note that regular stacks have real Z{0) by 
symmetry and hence P(0) — for any TV. We have 
ZxiO) ^ 1 at large positive T where neutral diagrams 
with g ~ dominate in Eq. JSJ and Zx{0) ~ — 1 at 
large negative F where ionic diagrams with q ~ ±7r dom- 
inate. Zx{0) vanishes at some F ~ 0, as seen for both 
iV = 14 and 16. The Jahn- Teller instability of An rings 
appears as a symmetry change of \['(0)S242i where Z{Q) 
vanishes and changes discontinuously in Fig. ^p. The 
arrow at Fc = —0.666 marks the NIT of the extended 
regular stack as found from extrapolations of symme- 
try crossovers. ^JS- Finite-size effects in Z{0) are reduced 
considerably in —N\n{\Z{Q)\'^), which is the proper size- 
independent quantityi^ The aei results in Fig. ^ nearly 
coincide for TV = 14 and 16 except around Z ^ 0. In 
Section 5 we report almost identical aei for = 14 and 
16 over the entire F range for stacks with equilibrium 
dimerization and Z ^ 0. 

Figure|21shows aei{6, F) for N — 16. Finite alternation, 
(5 > 0, strongly suppresses the divergence at the NIT, as 
seen on comparing S — 0.02 and (dashed line), and 
shifts the peak to F = at 5 = 1. Regular stacks with a 
finite potential against dimerization have finite polariz- 
ability for F > Fp > Fc. The polarizability of dimerized 
stacks with F < Fp remains finite because S > opens 
a gap in the energy spectrum. The S = 0.99 curve fol- 
lows closely the simple analytical result for S = 1 and 
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FIG. 2: Logarithm of the linear electronic susceptibility, aei, 
vs r for a 16 site stack with S =0.02, 0.05, 0.1, 0.4 and 0.99. 
The arrows mark the direction of increasing S, and the dashed 
line for S — — diverges at the NIT. 



independent dimers. 



aei = (r^ + 8)- 



(12) 



when there is no electronic delocalization. The large vari- 
ation of ttei {S, F) in Fig. El is understood in terms of lo- 
calization around F '--^ due to the dimerization gap and 
for |F| >> due to charge localization for any S. The 
pronounced asymmetry of a^i {S, F) with respect to F = 
is due to strong correlations in Hq(S,T) leading to quite 
different GS on the neutral and ionic sides. 

We illustrate in Fig. |2|the iterative solution of Eq. 1^. 
Circles refer to the first iteration with AM(0) for N = 16 
and 6 = 0.1. Crosses and squares refer to the second and 
third iteration, respectively. The results for aei = P'i^) 
do not change, since the linear perturbation in F appears 
from the outset and iterations introduce higher-order cor- 
rections. Finite S is required for finite P"{0) and cor- 
rections must be included. The second iteration more 
than doubles P"(0), but the third iteration produces no 
additional change. Taylor expansion of AM{F) through 
the linear term suffices for F^ contributions. Nonlinear 
polarizabilities based on the PBC formulation in Eq. (0) 
or raise interesting issues for rigid lattices. Here we 
focus instead on vibronic contributions in Eq. |j2Jl to the 
linear polarizability of models with Pcicrls transitions. 



IV. POLARIZABILITY NEAR THE PEIERLS 
TRANSITION 

The Su-Schrieffer-Heeger (SSH) models for polyacety- 
lene, (CH)^; , describes in the adiabatic approximation the 
Peierls instability of Hiickel chains with [/ = A = in 
Eq. H10|) and linear e-ph coupling Ue^ph — tS/u. Simi- 
lar approaches have been applied to Peierls transitions 
in segregated (A = 0) stacks of 7r-electron donors or 
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e,=0.28 



E =0.64 



FIG. 3: The linear electronic susceptibility, a^i, and the sec- 
ond derivative of the GS polarization on the electric field, P" , 
calculated as a function of F for a 16 site stack with 5 — 0.1. 
Circle, squares and crosses refer to results obtained at the first 
second and third iteration, respectively (see text). 



acceptorspi or to uncorrelated {U — 0) chains"^^ with site 
energies ±A. The stability of mixed stacks described by 
Hq has been studied through the response of a regular lat- 
tice to the SSH coupling. '^^ In the adiabatic approxima- 
tion, the GS potential energy surface is defined by adding 
the elastic energy for lattice motion to the electronic en- 
ergy. For harmonic potentials and linear electron-phonon 
coupling, the elastic energy of the Peierls mode is S'^ /2ed 
per site, where = o?f,_ph/k is the small polaron binding 
energy and k is the lattice force constant. ^■^ We therefore 
define the Peierls-Hubbard model by adding the lattice 
energy to Hq{5^T) in Eq. The electronic energy in 

an applied field is given by Eq. The total GS energy 
per site is 



8t^£o{5.T,F) 



2ed 



(13) 



t/,A 



The same development holds for arbitrary 
Eq. H1U|) or for other cell models with PBC. 

Eq. H13|) fully defines the GS potential energy surface 
with one additional parameter, the lattice stiffness l/e^, 
besides those entering the electronic Hamiltonian. The 
equilibrium dimerization is obtained by minimizing of the 
total energy, 



5eq — —fid 
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(14) 



The limiting inverse stiffness is = -\/2, which yields an 
equilibrium dimerization in Eq. (|14|l of (5 = 1 at F = 




FIG. 4: Equilibrium dimerization, S, softened frequency flp 
in units of the reference frequency cup, and infrared oscillator 
strength, calculated for Eq. (13) with — 0.28 and 0.64 
(left and right panels, respectively.) Circles refer to A'^ = 14, 
crosses to TV = 16. 



for a stack of decoupled dimers whose GS is easily found. 
The top panel of Fig. ^ shows the equilibrium dimeriza- 
tion of Ho{5, F) at F = for = 0.28 and 0.64 in units 
of t. Since results for N = 14 (crosses) and 16 (circles) 
now coincide almost exactly, we conclude that finite-size 
effects have become negligible in deformable lattices with 
S > 0.1 around F ~ that precludes strong dclocaliza- 
tion. 

Previous modeling of TTF-CA used the estimate Cd ^ 
0.15 - 0.3.2ii2^ The = 0.64 stack is less stiff, with 
a three-fold increase in the maximum dimerization to 
0. The Peierls instability occurs 
0, where the curvature l/e^ of the 
lattice potential and the curvature of the electronic en- 
ergy cancel exactly. The stack is regular for F > Tp{ed) 
and dimcrized for F < Fp(ed). As expected on general 
grounds, the neutral lattice with F >> is a band in- 
sulator that is conditionally stable. The divergence of 
{d^Eo/dS^ jo at Fc ensures a Peierls instability at some 
Fp > Fc that is model dependent. Dimerization de- 
creases but persists as F becomes more negativei^ The 
F << limit with p I maps into the spin-1/2 Heisen- 
berg antifcrromagnetic chain"^^ whose GS is unstable to 
a spin-Peierls transitioni^ 

The curvature of the GS potential energy surface is the 
frequency of the Peierls mode fip = {d^Sr /dS^)eq- At 
F = the ratio Qp/cup at the equilibrium (5 is a function 
of F: 



5eq ~ 0.3 around F 
at {d'^£T/dS^)o 



np 

UJp 



— = 1 + Q 
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(15) 



eq 



The middle panel of Fig. 2]shows the evolution of fip/wp 
for the correlated model HQ{6,r) in Eq. ((TT|l . The soft- 
ening of the Peierls mode on the neutral side and its sub- 
sequent hardening are observed for TTF-QBrCla whose 
continuous dimerization transition has been investigated 
by far-IR spectroscopyi^ii The temperature evolution of 
a combination band in the mid-IR spectrum of TTF-CA 
for T > Tc, also supports the presence of a soft lattice 
modet^ The data are in both cases consistent with uip ^ 
100 cm~^, a typical frequency for lattice vibrations in 
crystals with molecular masses of ^ 100 AMU. 

The Peierls mode is IR active by symmetry in mixed 
stacks and borrows huge IR intensity from electronic de- 
grees of freedomi^ In the present approximation, with M 
in Eq. Q for fixed sites, the IR intensity of the Peierls 
mode is entirely due to charge fluctuations induced by S. 
The corresponding oscillator strength is: 

= = —T- U J ^''^ 

Here me and e are the electronic mass and charge, respec- 
tively, |/i/_Rp is the squared transition dipole per site, and 
ed is in units oft. Once again the IR intensity is governed 
by P-derivatives and e^j. However ujp, a and t enter the 
expression for the dimensionless fip. The bottom panels 
of Fig. 0]show fjji for typical parameters^?^ a — 3.7 A, 
wp = 100 cm-i, t = 0.2 eV. 

The intensity of the Peierls mode has a pronounced 
peak at Fp where the electronic charges are maximally 
mobile. Dimerization localizes charges and lowers //^j 
for F < Fp. We have //p = when the electronic flux 
induced by dimerization reverses from D — > A on the 
neutral side to A ^ D on the ionic side. The neutral- 
ionic interface of dcformable lattices can be identified as 
fiR = 0, which occurs near F = Fc for regular stacks or 
Zx{S, F) = in dimerized stacksi^S The actual IR inten- 
sity does not vanish due to the motion of the molecular 
sites with charges ±ep, but frozen charge contributions^^ 
are neglected in Eq. Q . In materials that dimerize in the 
neutral side, the N-I interface could, in principle, appear 
experimentally as a dip in //p. However, this dip does 
not mark a phase transition. There is no change in sym- 
metry and all properties of the system vary continuously 
at the interface i^iiSS 

The IR intensity of the stiff (e^ = 0.28) stack in Fig. 
2] is an order of magnitude larger than that of the ed — 
0.64 stack. While contrary to the factor in Eq. (|15|l . 
this is readily understood in terms of greater delocaliza- 
tion in stiff lattices whose Fp approaches Fc and whose 
dP/dS diverges in the limit — ^ 0. At the Peierls tran- 
sition, the large oscillator strength of the dimerization 
mode corresponds to effective masses m* = rUe/ fiR ^ 
1000 and 6000 for ed = 0.28 and 0.64, respectively. So 
TO* at the transition is roughly a proton mass and is fully 
two orders of magnitude smaller than molecular masses. 
This justifies a posteriori the neglect of frozen charge 
contributions to //p near the Peierls transitioni^ 
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FIG. 5: Electronic and vibrational contributions to the static 
polarizability (dotted and continuous lines, respectively), for 
= 16 and the same parameters as in Fig. 0] The three 
lines for a^ib have, from top to bottom, for 'y/cop = 0.01, 0.1, 
0.2. 



The small effective mass associated with dimerization 
at the Peierls transition implies strong mixing of elec- 
tronic and vibrational degrees of freedom. Lattice con- 
tributions to the polarizability can be substantial, as an- 
ticipated for amb in Eq. To demonstrate this, we 
differentiate both sides of Eq. (|14|l with respect to F and 
use P — —{880/ dF) to obtain 

dSeq ^ ed dP_ 

OF 1 + edid^£o/dS^)eq as ^ ' 

We then use Eq. IjlSfl to write the vibrational polarizabil- 
ity of the Peierls mode at S^q' 




This expression is exact when the vibrational kinetic en- 
ergy is neglected, which is a good approximation in view 
of the low frequencies involvedi^S It is equivalent to the 
sum-over-states expressionism a^if, — 2\iijr\'^ jfiVlp, based 
on IR transition moments and energies. 

Large avib is expected for the dimerization mode due 
to its large transition dipole and low frequency. Indeed, 
57p — at the Peierls transition gives a divergent po- 
larizability. We suppress the divergence by introducing 
damping 7 that represents the lifetimes of lattice modes, 
anharmonic potentials, etc. Damping is introduced em- 
pirically by changing (fip/wp)^ as shown in the second 
equality in Eq. H18(l . Figure [S] reports a^b for the ed = 
0.28 and 0.64 systems shown in Fig. ^for damping ^jiop 
= 0.01, 0.1 and 0.2. The dotted line is the electronic 
polarizability P'{0) of the deformable lattice. The indi- 
cated damping corresponds to 7 = 1, 10 and 20 cm^^ for 
an estimated ujp = 100 cm^^. Quite predictably, damp- 
ing affects ayib in a narrow region about Fp and is more 
pronounced in the softer lattice. We note that typical 
bandwidths of lattice modes in molecular crystals are in 
the 1-10 cm"-'^ range. We will use 7 = O.lwp to model 
the dielectric response of CT salts and cannot specify 
ayib peaks more accurately than shown in Fig. |31 
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V. DIELECTRIC ANOMALY OF 
PEIERLS-HUBBARD MODELS 

We now combine the electronic and vibrational polar- 
izabilities of the Peierls-Hubbard model in Eq. to 
obtain the dielectric constant of the equilibrium lattice. 
In SI units, 



Here eo is the vacuum permittivity constant, Kqc '^-^ 3 is 
the usual contribution from molecular excited states that 
are not being modeled, and v is the volume per site. We 
are interested in remarkably large peaks with k > 100. 
To get proper dimensions for k the polarizabilities {aei 
and a.vib) are multiplied by e^a^/t. We adopt typical 
TTF-CA values, w = 206 A^, a = 3.7 A, and t = 0.2 eV, 
and in the lower panels of Fig. El we show the resulting 
K — 3 + 60(ae; + a„if,) for iV = 14 and 16 stacks with equi- 
librium dimerization and 6^= 0.28 and 0.64. The dotted 
lines are the electronic contribution to k. The dielectric 
peak is clearly associated with the Peierls transition and 
vibrations. The divergence of aei in Fig. ^at the NIT 
(Fc = -0.666) of the i5 = stack is strongly attenuated by 
dimerization and appears as a small bump around Fp in 
Fig. El The electronic fluxes leading to the IR intensity 
of the Peierls mode are also responsible for the large k 
peaks of 1000 and 100 estimated for = 0.28 and 0.64, 
respectively. We note that neglecting the softening, i.e. 
by imposing f2p — LUp in Eq. H18I) . lowers the k peak by 
an order of magnitude, while neglecting damping (7 = 0) 
gives a divergence. We took 7 = a;p/10 = 10 cm^^ in 
Fig. 6. 

Both the height and shape of the k peaks as a func- 
tion of temperature compare favorably to the dielectric 
data, discussed below, of Horiuchi et al&LSiSiiS That is 



FIG. 7: The equilibrium dimerization amplitude, S, the GS 
ionicity, p, and the dielectric response, k, for the same param- 
eters as in Fig. |l]except for V = 2. Circles and crosses refer 
to N =14 and 16, respectively, with continuous and dashed 
lines joining exact and mf results. Dotted lines in the bottom 
panels are the electronic mf dielectric response for A*' = 16. 

not the case for the ionicity p shown in the upper panel 
of Fig. |S1 We properly have p ~ 0.40 at the Peierls tran- 
sition of the ed — 0.28 stack and p ~ 0.25 in the softer 
stack, but p{T) hardly shows any sign of a transition. 
By contrast, the measured ionicity has a kink in systems 
with continuous transitions around p ^ 0.25&2i2424iS The 
basic model Ho{S,r) in Eq. Hll|l is deficient in this re- 
spect, just as its continuous NIT in rigid lattices can- 
not account for discontinuous p at Tc ^ 81 K in TTF- 
CA. Several extensions of the model have long been rec- 
ognized. Coulomb interactions occur both within and 
between stacksi224£ Coupling to molecular modes is an- 
other general phenomenoui^' In practice, mean-field (mf ) 
approximations are necessary for inter-stack interactions. 
While observing a p anomaly at the Peierls transition re- 
quires going beyond Hq (S, F) , the appropriate extensions 
and parameters remain open. 

To illustrate a simple extension that is suitable for ei- 
ther exact or mf analysis, we introduce nearest-neighbor 
Coulomb interactions V along the stack, 

Hv = Hoi6,r) + VJ2<lp<lp+i+N^ (20) 

Since V is diagonal in the VB basis, all GS properties 
arc found as before. The solid lines in Fig. [7|are exact 
A'^ = 14 and 16 results for the equilibrium dimerization, 
ionicity and dielectric anomaly of stacks with V — 2 and 
otherwise the same parameters as in Figs. 01 and El The 
dashed lines are the mf approximation to the V term in 
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FIG. 8: The same as Fig. 13 for V = 3; for the sake of clarity 
only mf results are reported, with circles and crosses referring 
to N = 14, and 16, respectively. Dotted lines in the lowest 
panels show the electronic contribution to «. 



Eq. (|2()|l . The mf results for V — 2 are similar, with 
somewhat sharper p at Tp as expected for an interac- 
tion that stabilizes adjacent ion-pairs. Both = 0.28 
and 0.64 now have a kink at p{Tp) that qualitatively re- 
sembles p{T) data. The k peak is lower by a factor of 
two, while the maximum dimerization hardly changes. 
Charges are kept together for V > 0, thereby reducing 
the polarizability. 

The first-order transition of TTF-CA at 81 K has 
a Ap 0.2 and concomitant neutral-ionic and Peierls 
transitionsiSi^ The ionicity is a suitable order parame- 
ter in this case. Discontinuous p has long been treated in 
mf theory^SS and mf results for V~ 3 are shown in Fig. 
|S1 for the same quantities and parameters as in Fig. |7| 
Now ed = 0.28 produces a Ap jump in stacks of A'^ = 
14 and 16 that resemble previous results for smaller N 
up to 12i^ The kink in p for ed = 0.64 is stronger and 
still larger V will make it discontinuous. The maximum 
dimerization remains almost independent of V. The k 
peak now appears at the NIT and is reduced for both 
ed — 0.28 and 0.64. The reason for reduced vibrational 
contributions is that a discontinuous NIT interrupts the 
softening of wp. This is shown in Fig. |^for the same 
parameters, together with //p. 

Figure [3 contrasts exact and mf results for small V = 
2, while Fig. |S1 shows a first-order mf transition for V — 
3 in the stiff stack. To summarize the various roles of V 
in Eq. (|20|l . we note that in mf it also accounts for the 
coupling of electrons with molecular vibrations. Totally- 
symmetric (ts) molecular vibrations modulate on-site en- 
ergies, and hence F, leading to Holstein coupling^i The 



effects of this coupling on vibrational spectra of CT salts 
are well knowniiSiiSi^ In regular stacks, ts vibrations 
only appear in Raman spectra, whereas in dimerized 
stacks they acquire large IR intensity, proportional to 
[dP/dVY, through their coupling to electronic degrees 
of freedom4^ In the adiabatic approximation, Holstein 
coupling plays exactly the same role as V within mfn^ii^ 
Thus mf results in Figs. E|and|Hlrepresent an effective pa- 
rameter V whose interpretation is model dependent. Al- 
though the large IR intensity of Holstein modes in dimer- 
ized lattices suggests that also they contribute to k, their 
high frequency 1000 cm"-'^) makes this contribution 
negligible in comparison to the Peierls modes. 

VI. DISCUSSION 

The dipole P for insulators with PEC and its ex- 
tension in Eq. ^ to an applied field, make it possible 
to model the electronic and vibrational polarizability of 
quantum cell models. We found the induced-dipole oper- 
ator AM in Section 2 for quantum cell models in general 
and applied it to a modified Hubbard model, Hq{5^T) 
in Eq. 1)11(1 . We report model-exact electronic polariz- 
abihties in Fig. |21 for the NIT of regular {5 = 0) and 
dimerized stacks. We then include the Peierls mode with 
frequency wp in the adiabatic approximation and its soft- 
ening. The IR intensity of the Peierls mode in Fig. 0| 
and the dielectric peaks in Figs. El and [T] are related to 
Peierls transitions in stacks with stiffness 1/ed- We have 
chosen — 0.28 as a plausible estimate for TTF-CA and 
a softer stack with ed= 0.64 whose Peierls transition is 
around p ~ 0.3 on the neutral side. The choice of < = 
0.2 eV comes from TTF-CA optical data. The model pa- 
rameters are not arbitrary, but typical of values used for 
other propertiesj^i*^ The same analysis holds for various 
extensions of Ho{S,r) for organic CT salts. All results 
are derived from the F-dependence of the GS properties 
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of stacks described by Ho{5,T), and at the equilibrium 
dimerization for fixed in Eq. and fixed nearest- 
neighbor interactions V in Eq. (|2()(l . 

The dielectric data of Horiuchi et al. are taken as a 
function of temperature at ambient pressure ^iSiSiiS or, 
more recently, as a function of pressure, up to ~ 10 kbar, 
at constant temperatureiii Decreasing volume on cool- 
ing or compression stabilizes the ionic phase by increas- 
ing the crystal's Madelung energy. Cooling corresponds 
to decreasing F and increasing p in models, but the rela- 
tion between temperature or pressure and F is not known 
at present. Cooling also increases DA overlaps that en- 
ter in t and for one-dimensional models. Overlap be- 
tween stacks also increases, but although the importance 
of inter-stack interactions has been notedj^ their model- 
ing is still rather qualitative. 

The crystals studied by Horiuchi et al. in refs. 0, O 
have D = DMTTF (dimethyl-TTF) and the quinones A 
= QBr„Cl4_„ with n = 0-4; n = and 4 are CA and 
bromanil (BA), respectively. The acceptors have two- 
fold disorder in Br/Cl in the crystal lattice for n = 1, 2 
and 3. There are two QBr2Cl2 isomers, a centrosymmet- 
ric one called 2,5 and a polar one called 2,6. The n = 1, 3 
and 2,6 complexes are formally polar, but only weakly so: 
Their dielectric peaks are similar but are shifted to lower 
temperature. As there is no disorder in the prototypical 
system, TTF-CA, the best candidates with continuous 
transitions are DMTTF-CA and DMTTF-2,5QBr2Cl2. 
The dielectric peak of DMTTF-CA at 65 K is k - 220 
at 30 kHz; p increases quickly from ~ 0.32 to ^ 0.42 
between 65 and 60 K and then slowly to ^ 0.48 at 10 
K. The DMTTF-2,5QBr2Cl2 peak at 50 K is k - 170 at 
30 kHz, with p increasing from 0.29 to 0.36 on cooling 
from 50 to 10 K. The n peak of DMTTF-QBrCls at - 55 
K closely resembles the 2,5 system, but p has not been 
reported. Pressure-induced k peaks are reported^-'^ at 5 
K in DMTTF-BA {n = 4) and in DMTTF-2,6QBr2Cl2, 
both with p < 0.3 at 5 K and 1 atm. 

The observed peak heights of k are consistent with di- 
electric peaks of the model at the Peierls transition. We 
conclude that a n peak marks the Peierls transition in 
systems with continuous p. The narrower k peak in Fig. 
13 for model parameters leading to a first-order transition 
is also seen in the narrow k — 600 peak of TTF-CA. * 
The peak shapes and heights follow qualitatively the cal- 
culated pattern of increased polarizability in stiff lattices 
with Peierls transitions (continuous p) occurring at Fp 
close to Fc. More quantitative analysis requires a relation 
between temperature and F. As illustrated in Section 
5, accurate modeling of p{T) near Fp depends on such 
specifics as damping and intra and inter-stack interac- 
tions. By contrast, the Peierls transition of deformable 
stacks and NIT of rigid stacks are general features of 
Hubbard and related models. 

The dielectric anomaly observed in mixed stack CT 
salts is clearly associated with lattice degrees of free- 
dom and the Peierls transition in systems with either 
continuous or discontinuous NIT. The huge IR intensity 



of the Peierls mode around the Peierls transition^^ and 
its softening are responsible for ayib >> ctei in Eq. 
Available spectroscopic dataiS in the far-IR region for 
TTF-QBrCls support this picture.^^ It is interesting to 
compare the Peierls transition occurring in this system 
at p ~ 0.3 with the discontinuous NIT occurring in TTF- 
CA at similar p: in TTF-CA the NIT occurs before the 
complete softening of the Peierls mode and a reduced Hi- 
peak is observed. Whereas far-IR data are not available 
for TTF-CA, the incipient softening of the dimerization 
mode has been recently extracted from a detailed study 
of combination bands in the mid-IR region.''^ A more 
systematic analysis of far-IR spectra and dielectric prop- 
erties of mixed stack materials is certainly desirable as 
it will confirm the important connection between lattice 
degrees of freedom and materials properties. 

Organic CT crystals have DA repeat units along the 
stack and, in contrast to most inorganic salts, are quasi- 
one-dimensional systems. Rapid convergence, often as 
is typical for GS properties of systems with PBC 
and small repeat units. It is very convenient for model- 
ing that P{F) is a GS property. As found throughout, 
all = 14 and = 16 results coincide in deformable 
lattices with 5 > 0.1 at F ~ 0. Larger N is not required 
for CT salts, but is needed for the NIT of (5 = stacks. 

Discontinuous p at the NIT is fundamentally differ- 
ent from continuous p at the Peierls transition of softer 
lattices. We have focused on the Peierls-Hubbard model, 
Hq{5, F) in Eq. Hll|l . which has continuous p and a Peierls 
instability on the neutral side for large e^. The order 
parameter S for the second-order transition breaks the 
rcfiection symmetry of the regular chain. Once the stack 
has dimerized, there is no further transition and, indeed, 
hardly any remnants of the prominent neutral-ionic tran- 
sition of the 5 = stack. Models with continuous p have 
Peierls rather than neutral- ionic transitions. Similar re- 
sults are obtained for stacks with not too large y-like 
interactions: The Peierls or structural transition occurs 
at small p, and V promotes a more prominent kink in p. 
Increasing further V (or Holstein) interactions leads to a 
discontinuous NIT, where the on-site charge reorganiza- 
tion is the driving event, p is the proper order parameter, 
and dimerization just follows from the unconditional in- 
stability of the ionic phase. 

The increasing number of CT salts undergoing phase 
transitions demonstrates^^ far more diversity than the 
neutral-regular and ionic-dimerized GS of TTF-CA. The 
lattice stiffness l/sd is an important parameter for dis- 
tinguishing between Peierls and neutral-ionic transitions. 
Disordered D or A in single crystals raises different is- 
sues. The microscopic modeling of CT salts has primarily 
been in terms of Hubbard-type models with many open 
questions about parameters. We anticipate that accu- 
rate modeling of P{F) and its applications to dielectric, 
vibrational and other data will lead to more quantita- 
tive descriptions of quantum phase transitions in these 
materials. 
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